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Abstract 

We study the geometry of propagation of relativistic fermions. We propose how to measure 
its quantum mechanical length. Numerical lattice results for the free propagator of Dirac- 
Wilson fermions yield Hausdorff dimension rfj^ = 2 for the unit-matrix component and 
dn = ^ for any 7-matrix component. A possible generalization when matter interacts with 
radiation is discussed. 
PACS index: 03.65.-w, 05.30.-d 







1. Introduction 

What do we know about the geometry of propagation of massive particles in non-relativistic 
quantum mechanics? According to Feynman and Hibbs quantum mechanical paths are 
zig-zag lines, which are no-where different iable, exhibiting self-similarity when viewed at 
different length scales. In terms of modern language this is a fractal curve [0]. Abbot 
and Wise have shown for free motion that quantum mechanical paths are fractal curves 
of Hausdorff dimension dn = 2. Actually, in quantum mechanics the concept of paths is 
not well defined. If, however, one goes over from real to imaginary time the path integral 
becomes mathematically well defined (Wiener measure). In imaginary time (Euclidean) 
quantum mechanics the free motion resembles the Brownian motion of a classical particle 
Then one has a stochastic interpretation as classical Brownian motion and paths can 
be considered as random variables of a Gaussian process The typical path of a classical 
particle carrying out a Brownian motion is a fractal curve with dn = 2. According to 
Itzykson and DroufFe when following a Brownian curve on a lattice from xq, to to Xi, ti 
the typical distance behaves as 

|xi-fo| ~(ti-to)", ^^ = 1/2, (1) 

and u is called a critical exponent, (by analogy with the power laws in the theory of critical 
phenomena). Considering the typical distance as a function of the elementary length Ax 
then defines the Hausdorff dimension (see Eq.(2) below). It is important to note that 
dn = 2 for free motion in quantum mechanics does not change when going over from real 
time to imaginary time. 

For the case of a massive interacting quantum mechanical particle Campesino-Romeo 
et al. have shown analytically dn = 2 for paths of the harmonic oscillator. Other 
potentials have been investigated by numerical simulations on the lattice pj. It turned 
out that dn = 2 for local potentials (harmonic oscillator. Coulomb potential). However, 
dn 7^ 2 was found for velocity-dependent interactions. Such velocity-dependent actions 
are supposed to play a role in condensed matter (dispersion relation of a massive particle 
travelling in a medium) or in Brueckner's theory of nuclear matter |p. 
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The Hausdorff dimension of field configurations has recently become a subject of interest 
in quantum gravity (in particular in 2-D), where the role of fields is played by the geometry 
of space-time - [0. The Hausdorff dimension can, e.g., be defined as a power law 



relation between two dimensionful observables at a critical point [|T0|- E-g-, for quantum 
gravity in D=2, one has dn = 4: [0.1t is interesting to establish relations between the 
Hausdorff dimension and critical exponents. Such a relation is known to exist between dn 
and the critical exponent z/, which is the analogue of a critical exponent of the spin-spin 



correlation function in statistical mechanics. In quantum gravity |T6| z/ is defined as the 
critical exponent of a mass, m(A/i) ~ {AfiY, when the cosmological constant tends to 
its critical value A/i 0, where the mass m(Ayu) characterizes the fall-off behavior of a 
two-point function G^{r) ~ exp[— m(Ayu)r]. Then holds the following relation (scaling law) 



dn = l/i^, relating the Hausdorff dimension to the critical exponent v [|^. Thus also du 
plays the role of a critical exponent. 

Critical exponents of a theory determine its universality class and thus classify the the- 
ory. The above example of quantum gravity leads us to the question: Can the Hausdorff di- 
mension play a similar role for models like, e.g., theory of matter interacting with radiation, 
in condensed matter, in the theory of medium energy nuclear physics (Dirac phenomenol- 
ogy), or in high energy for the theory of quarks and gluons (plasma)? A related point of 
view is the following: In quantum mechanics the Hausdorff dimension of typical paths may 
differ from the standard value d// = 2 in the presence of velocity-dependent potentials 
also if one considers quantum mechanics in a background medium corresponding to curved 
space-time. Such a situation occurs, e.g., when considering the relativistic propagation of a 
particle in nuclear matter (e.g., a neutron star), or when considering a relativistic particle 
impinging on the surface and propagating a short distance in ordinary matter. 

With those questions in mind we want to start out here by asking the much simpler 
question: What happens to the geometry of propagation of a massive particle in relativistic 
quantum mechanics? We want to study the question by considering the most simple massive 
relativistic particle occuring in nature: the non-interacting Dirac fermion. The purpose of 
our work is firstly to compute du in general and secondly to suggest a definition based on 
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the propagator which allows for a generalization to interacting fermions. 



2. How to measure the geometry of paths in non-relativistic quantum 
mechanics? 

The Hausdorff dimension of a fractal curve is defined as follows: Suppose one has an elemen- 
tary length scale (resolution) Ax to cover the curve. Experimentally, Ax corresponds to 
the resolution of an experimental apparatus (yardstick, wavelength of light in microscope). 
Measuring the length L of a curve in terms of an elementary length Ax, the property of 
being fractal is captured in the Hausdorff dimension d^, defined by 

Lo~,^o^e'""\e = Ax/Li (2) 

where Li is a fixed length. The Hausdorff dimension dn is a number chosen such that Lq 
becomes independent of e in the limit e ^ 0. Hausdorff has given a precise definition of 
"resolution" by covering the cutrve with L/Ax spheres of diameter Ax. 

The definition given by Eq. ( g) has been applied in Ref . P, || to characterize the typical 
path of a quantum mechanical trajectory and in this sense the Hausdorff dimension has 
been computed. But in order to be precise one has to say what is meant by length L and 
resolution Ax for quantum trajectories. Suppose we consider the amplitude for propagation 
of a particle, being at t = in a state characterized by position Xin and at t = T in a state 
characterized by position x/j. We discretize time to = < ti < ■ ■ ■ < = T, with At = S 
and T = NAt. Then the amplitude is given by a path integral, which after discretization 
of time reads 



Z{5) = J dxi---dx]\! exp[iS[xk, S]/h] 



(3) 

We have denoted x^ = x{tk). This amplitude is an approximation of the continuum (exact) 
amplitude, obtained by taking the limit At 0, i.e., ^ oo in an appropriate way. A 
suitable observable to study the geometry is the propagator length defined by [§] 

< L{S) > = <^|xfc+i-Xfc|> 

fc=0 



Xq — X'lji ,X ]\j — fi 

■ (4) 

XQ — X'Iyi^^N — fi 



1 f 

— / dxi---dxN \xk+i-Xk\ exp[iS[xk,S]/h] 



In quantum mechanics, the resolution of length is given by the dynamics. It is natural to 
define it as average length increment corresponding to a time increment 6 0, 

N-l 



<A.>^i<i:i..«-..i> 



(5) 

XQ — X^Yi j\f — fi 



As fixed length Li one can choose, e.g., the length of the classical trajectory from x^, t = 
to Xfi,t = T given by the classical continuum action. If one is interested only in dn, i-e., 
the exponent of the power law Eq.( 0), then any fixed length Li can be chosen, in particular 
the length of the straight line between Xm and x/i. In the continuum limit of quantum 
mechanics one takes the limit At = 6 —>■ 0, which implies < Aa; >~ e ^ 0. The length 
< L > has been computed by numerical simulations on the lattice using imaginary time 
(Euchdean) quantum mechanics and dn has been extracted M. 



3. Definition of propagation length in relativistic quantum mechanics 

When trying to generalize the above considerations to relativistic quantum mechanics, one 
is faced with the following problem: (a) In quantum mechanics, position of a particle is an 
observable and can be measured. In quantum field theory position is not an observable, 
(b) The particle number is conserved in quantum mechanics. In quantum field theory 
particles can be created and annihilated, (c) In quantum mechanics particles propagate only 
forward in time. In quantum field theory particles propagate forward and backward (anti- 
particles) in time. Thus in a relativistic theory one must look at space and time dependence, 
corresponding to a causal propagation of a massive particle. From the mathematical point of 
view as we have pointed out above, the path integral in imaginary time quantum mechanics 
is well defined in terms of a stochastic process. This can be generalized to Euclidean 
path integrals of bosonic (polynomial) field theory which are mathematically well defined, 
allowing an interpretation as stochastic process H. The measure gives the dominant 



contributions of no- where differentiable curves [18|. However, for fermion field theory, there 



is not yet a strict mathematical formulation in terms of a stochastic process, although 



Osterwalder and Schrader 19| have established a Feynman-Kac formula for fermion fields. 
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From those remarks it is evident that the definition of propagator length given by Eq. (4) 
for non-relativistic quantum mechanics can not be simply taken over to relativistic quantum 
mechanics. The definition (4) is based on a discretization of time (taking At = 5 finite and 
letting 5 ^ in the end when extracting the critical exponent). In relativistic quantum 
mechanics we discretize time and space, i.e., we work on a lattice with some finite lattice 
spacing a (letting a — in the end). Let us introduce a new definition for the length of 
propagation for relativistic quantum mechanics, by considering, in particular, the Dirac 
fermion action. It is well known that the corresponding lattice action is plagued by the 
so-called fermion species doubling problem: 16 copies of fermions with the same mass occur 
(as poles of the fermion propagator). A way out commonly used is the Wilson action, which 
lifts the species doubling by adding another term. The (Euclidean) Wilson-Dirac action on 
the lattice reads [^] 

S[i),^] = '^i^mKm,n'ipn- (6) 

m,n 

Here m,n are indices (tupel) which denote lattice sites (e.g., x„ = na, n = 0, ±1, . . . , in 
D = 1). The fields have been rescaled a^^'^(am + ArY^'^ip ^ ip, k = l/(2ma + 8r) is the 
hopping parameter and r is the Wilson parameter (usually chosen to be r = 1). The matrix 
K is expressed in terms of the hopping matrix M 
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Mm,n = '^{r + -ff,)Sm+/l,n + {r - J^)Sm-fi,n- (7) 

Here m + fi denotes the lattice site next to site m in the positive //-direction. The fermion 
propagator is given by 

< Ipni^rn >= {K-^)n,m- (8) 

The matrix M, coming from the kinetic term of the fermion action and from the Wilson 
term, allows the fermion to hop from one lattice site to the next neighbour lattice site. We 
suggest to define the length of the fermion propagator by counting in a non-perturbative 
way the number of hoppings. In particular we suggest the following definition 



The indices m, n denote the lattice sites, where the fermion is created and annihilated, 
respectively. By expanding the right-hand side of Eq.(9) as a power series in the hopping 
parameter, one finds that the p-th power of the hopping matrix M, which allows the fermion 
to hop a distance pa, gets multiplied with a factor p. Thus < > can be interpreted as 
a (non-perturbative) counter of how many times a fermion hops between sites m and n. 

The classical length is Ldass = \xm — Xn\- One has to compute numerically on the lattice 
< L > as a function of a/ Ldass- The goal is to extract a critical exponent 7, 

<L> 



[a/Ldass)-^. (10) 



Ldass 

By Eq.(2), 7 is related to dn via = 1 + 7. Because action (6) is parametrized in terms 
of the hopping parameter k, it is natural to consider another critical exponent a defined 

by 



^ ^ ^K^KcHt ( ^^^^ J ■ (11) 

^class ^ i^crit ' 

The critical exponent 7 is defined in the continuum limit a — > 0. When a goes to zero, the 
dimensionless lattice mass ma goes to zero, and k goes to its critical value Kcr-it — ^/'^t. For 
a free Euclidean fermion theory, which we investigate numerically, both exponents coincide, 
q; = 7. 

4. Numerical results 

Before measuring the length of the fermion propagator, it is useful to see if the definition 
of the propagator length makes sense. In order to have a meaningful length definition, one 
would expect < L > to obey a power law (11) when approaching Kcrit for fixed Ldass- In 
order to get a first idea on the behavior of < L > we have considered a drastically simplified 
hopping matrix Mj ,,• = + ^i+ij- We have dropped any dependence from 7-matrices 
and study the length as a function of k. We have done numerical calculations for D = 1 
with Kcrit = 1/2 on lattices up to = 120 and for D = 2 with Kcrit = 1/4 on lattices up 
to TV = 50. The numerical results confirm the expected scaling behavior of Eq.(ll). the 
results yield a — 0.49 for both D — 1 and D — 2. 
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Now we turn to the fermion propagator. For the free fermion case, Kctu = l/8r in 
D = A. The Euchdean free fermion propagator for Wilson fermions is given in momentum 
space by 

Afc = 1 - 2rK ^ COS A;^ + 2fi: ^ sm fc^ . (12) 
This is related to the space-time propagator A^^y =< ipx'^y > by Fourier transformation 

A... = :^Ee^'-^^"'^^,, (13) 
^ k 

where V = NiN2N^N4 is the lattice volume. Note that has a pole at A; = 0, Hcru = 
l/2Dr in D space-time dimensions. We choose periodic or anti-periodic boundary con- 
ditions. They correspond to the following choice of lattice momenta (see Ref. |pO||), 
kfj, = 2'7in^/Nfj_ corresponds to periodic boundary conditions and fc^ = 27r(n^ + 1/2) /N^ 
corresponds to anti-periodic boundary conditions. In both cases, G 0, 1, ■ ■ ■ , A^"^ — 1. 
Thus in the anti-periodic case /c^*" = /N^ is the minimal value of k^. 

We consider two components of the propagator: the unit-matrix component is given 
by \Tr[/S.x,y] and the 7^-component by \Tr[^^/S.x,y\- Let us consider the cases of D = 1, 
D = 2 and D = A space-time dimensions. For D = 1, one can compute the large lattice 
limit {V 00, not the continuum limit) analytically and obtains 

-^Tr[/\x—y=n\ = 2 K , 1 n, K ^ f^criti 

< L>= Lclass- (14) 



However, we are interested in the continuum limit a ^ 0, which corresponds to k — >■ 
Kcrit, with Kcrit = l/2r ill D = 1. The Wilson parameter is r = 1. The numerical 
results for the 71-component and the unit-component are shown in Fig. [la,b]. For the 71- 
component, we have chosen anti-periodic boundary conditions. We have varied N = Ni = 
4,8, 16, • ■ • , 1024. In order to approach Kcrit we have varied k = l/[2rcos(A;) -|- 2sin(/c)], 
with k = A;™*" = tt/N. We have chosen as classical length Ldass = \x — y\ = 1. We 
have evaluated the space-time propagator (13) by doing the Fourier transformation of A^ 
and of ^Afc numerically. From that we have evaluated the length < L > via Eq.(9) and 



hence the exponent a and via Eq.(ll). The resuh shows a power law behavior (11) with 
a = —0.0016, which corresponds to = 0.9984, by Eq.(2). 

The behavior of the unit- component, with periodic boundary conditions, is different. 
We have varied = 20, 40, • • • , 100. Because for periodic boundary conditions /c"**" — 0, 
we have approached k — > Kcrit by choosing k — 0.45, 0.475, 0.4875, • • • (decreasing \Kcrit — >A 
by a factor 2 in each step.). Now we have considered the classical length L^ass — N/2. The 
computation of < L > and dn is as above. One observes (Fig. [lb]) a power law with the 
critical exponent varying between a — 1.000 and a — 0.9983, and the corresponding fractal 
dimension varying between dn — 2.000 and dn — 1.9983. Generally, one observes that the 
larger the size of the lattice, the closer one has to be at Kcrit before the scaling behavior 
(11) is seen. How can a curve in topological dimension D = 1 show a fractal dimension 
larger than 1? The Hausdorff dimension measures the hopping of the fermion forward and 
backward on a line (D = 1), which can be fractal. 

Similar results are obtained m D = 2 space-time dimensions, shown in Fig.[2a,b]. Now 
Kcrit = l/4:r. Firstly, we have considered the 72-component. Because we have a regular, 
symmetric lattice, the ^-dependence is the same for all 7^-components. Thus we can 
interpret the 71-component as space component and the 72-component as time-component. 
We have chosen boundary conditions periodic in space and anti-periodic in time. We have 
varied A^i = Ngpace — 4, 8, 16 independently from N2 — Nfime = 4, 8, 16, • • • , 1024. In order 
to approach Kcvu we have varied k — l/[2r(l-|-cos(A;))-|-2 sin(/c)], where k — A;^*" = T^/Nume- 
As classical length we have chosen Ldass — Ntime/'^- We have obtained a — 0.0022 and 
dff = 1.0022. 

For the unit-component, we have chosen periodic boundary conditions in space and 
time. We have varied = A^i = A^2 = 10, 20, ■ ■ ■ , 50. In order to approach Kcrit we have 
chosen k = 0.225, 0.2375, ■ ■ ■ (decreasing \Kcrit — by a factor 2 in each step). As classical 
length we have chosen L^ass = N/2. We find dn = 1.999 to dh = 1.994 for lattices varying 
between = 10, • • • , 50. 

Finally, we present the results for D = 4 in Fig.[3a,b]. Now Kcrit = l/8r. For the 74- 
component, we have chosen boundary conditions periodic in space and anti-periodic in time. 
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Wc have varied A^^ = A^2 = A^3 = Nspace = 4, 8, 16 and = Nume = 4, 8, ■ ■ ■ , 256. In order 
to approach Kent we have varied k, = l/[2r(3+cos(/i;))+2 sin(/i;)], where k = k^^"" = Tr/Nume- 
As classical length we have chosen Ldass = Ntime/'^- As results we obtain a = 0.0086 and 
dn = 1.008. 

For the unit-component, we have chosen periodic boundary conditions in space and 
time. We have varied A?" = A^i = • • • = A^4 = 4, 8, 12 • • • , 28. In order to approach Kcru 
we have chosen k — 0.100, 0.110, • • • (decreasing \Kcrit — k| by a factor 2 in each step). As 
classical length we have chosen Ldass — N/2. We find dn — 2.10 to dh — 2.05 for lattices 
varying between A?" = 4, • • • , 28. 

One obtains the following picture: The 7^-conmponent of the propagator shows results 
compatible with dn — 1, i.e., no fractal behavior, in D = 1, 2, 4 and different combinations 
of periodic/anti-periodic boundary conditions. However, the unit-component shows results 
compatible with du = 2 for = 1, 2, 4, i.e., the same fractal behavior as in non-relativistic 
quantum mechanics. The numerical results, that is rf// = 1 for the unit-component and 
dn = 2 for 7^-component are independent of these boundary conditions. For larger lat- 
tices, scaling sets in later (closer to Kcrit)- Numerical errors increase when approaching the 
singularity k — > Kcrit- Also, numerical errors increase with the size of the lattice. Neverthe- 
less, one observes for D = 4 that when increasing the lattice size, the Hausdorff dimension 
moves closer to the value 2. 

5. Discussion 

(a) Why differ the results of dn for different components? Let us consider periodic bound- 
ary conditions and compare the unit-component with the 71-component of the propagator. 
Then the propagator, given in momentum space by Eq.(12), projects under the Fourier 
transformation (13) onto the cos(A; ■ (x — y)) part for the unit-component, but onto the 
sin(A; • [x — y)) part for the 71-component. In other words, the unit- and the 71-component 
have a different pole structure when k ^ Let us consider the continuum limit, a — >^ and 
e = Kcrit— >^ ~^ 0, but keep the lattice volume V fixed. Also we keep LcUss — \x~y\ — const. 
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and consider in D = 1 the unit-component ^Tr[Ax^y\, given by Eqs. (12,13), as a func- 
tion of e. For small enough e, the dominant contributions come from lattice momenta 
ki with ki << e. Taking in Eq.(13) only those lattice momenta into account, one finds 
lTr[A^J ~ 1/e. Consequently, K£\Tr[A^J ~ l/e^. Thus < L >~ 1/e, which by 
Eq.(ll) implies dn = 2. Doing the analogous calculation for the 71-component yields 
^Trl^iA^^y] ~ const, and n-j^^TrljiA^^y] ~ const.. Thus < L >~ const, and hence 
dn = ^- This is in agreement with our numerical results. 

(b) The results for the free fermion propagator on the lattice can be compared with 
Feynman's analytical expression for the fermion propagator in the asymptotic regime 
« and » Feynman expresses the propagator kernel i^+(2,l) = 



iii'Jfj.d^ + m)I+(t,x), and gives for the function /+ the asymptotic expression 

x) exp{-im[t - {x'^/2t)]}, x^ « (15) 

It can be seen that the propagation kernel is essentially the same as for a free particle in 
non-relativistic quantum mechanics, where the Hausdorff dimension is dn = '2- This is in 
accord with our result dn = "2 for the unit-component of the fermion propagator, which 
dominates the non-relativistic regime. 

(c) The definition of length (9) for the action (6) can be generalized to the case when 
matter interacts with radiation, i.e., QED. Then the fermion-photon interaction has the 
same structure as in Eq.(6), but the hopping matrix M[U] depends now on the gauge 
field via the link variables U^{n) (for details see [^), and the matrix element {K~^)m,n 
occuring in the fermion propagator, Eq.(8), must be replaced by a quantum expectation 
value < {K[U]'^)m,n > which means doing a path integral over the gauge field. However, 
one has to fix the gauge (see p2|). 

In summary, we have suggested a definition of length for the propagation of relativistic 
fermions. It shows scaling behavior when approaching the continuum limit and yields 
the critical exponents dfj = 2 (unit-component) and du = 1 (7^-component). This is 
consistent with the analytical behavior of the Fermion propagator. Our length definition 
can be directly generalized to interacting theories, e.g., matter with radiation. 
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Figure Caption 

Fig.l < L > /Lciass versus {Kcrit — i^)/K'crit for free fermion propagator in D=l, (a) unit 
component, (b) 71-component of propagator. 

Fig.2 Same as Fig.[l] in D=2 dimensions, (a) unit-component, (b) 72-component. 

Fig.3 Same as Fig.[l] in D=4 dimensions, (a) unit-component, (b) 74-component. 



14 



0.3 



0.2 



i 0.1 



V 



0.0 



■0.1 



-0.2 



1.0 



2.0 



D = l 



• djj = 0.998 



3.0 4.0 
- ln[(K - K)/K ] 



5.0 



6.0 




- ln[(K^ - K)/Kj 



0.7 - 



0.5 



0.3 



0.1 



-0.1 




1.0 



2.0 



3.0 4.0 5.0 

- ln[(K - K)/K ] 



6.0 



7.0 




- ln[(K^ - K)/Kj 



1.5 



ctl.O 



A 
V 



0.5 



0.0 




1.0 



2.0 



3.0 4.0 5.0 

- ln[(K - K)/K ] 



6.0 



7.0 




- ln[(K^ - K)/Kj 



